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Abstract 

We compute the supersymmetric contribution to Rb = T(Z — > bb)/T(Z — > 
hadrons) in a variety of supersymmetric models. In the context of supergravity 
models with universal soft-supersymmetry-breaking and radiative electroweak 
breaking we find i?^ usy ^ 0.0004, which does not shift significantly the Standard 
Model prediction (Rf M = 0.2157 for m t = 175 GeV; i?f p = 0.2204 ± 0.0020). 
We also compute Rb in the minimal supersymmetric standard model (MSSM), 
and delineate the region of parameter space which yields interestingly large val- 
ues of Rb- This region entails light charginos and top-squarks, but is strongly 
restricted by the combined constraints from B(b — > sj) and a not-too-large 
invisible top-quark branching ratio: only a few percent of the points with 
R s b usy > 0.0020 (la) are allowed. 
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1 Introduction 



Precision tests of the electroweak interactions at LEP have provided the most sensitive 
checks of the Standard Model of particle physics. The pattern that has emerged is 
that of consistent agreement with the Standard Model predictions. This pattern 
seems to have so far only one apparently dissonant note, namely in the measurement 
of the ratio Rb = Y(Z — > bb)/T(Z — > hadrons), where the latest global fit to the LEP 
data (-R^ XP = 0.2204 ± 0.0020 |3J) lies more than two standard deviations above the 
one-loop Standard Model prediction [fj] for all preferred values of the top-quark mass 
(e.g., Rf M = 0.2157 for m t = 175 GeV). Further experimental statistics will reveal 
whether this is indeed a breakdown of the Standard Model. In the meantime, it is 
important to explore what new contributions to Rb are expected in models of new 
physics, such as supersymmetry. 

The study of supersymmetric contributions to Y(Z — > bb) has proceeded in 
two phases. Originally the quantity e b was defined as an extension of the e 12 ,3 
scheme 0] for model-independent fits to the electroweak data. More recently it has 
become apparent that the ratio Rb is more directly calculable |5|, ||, 0] and readily 
measurable. It has been made apparent || that supersymmetric contributions to Rb 
are not likely to increase the total predicted value for Rb in any significant manner, as 
long as typical assumptions about unified supergravity models are made. On the other 
hand, if these assumptions are relaxed, it is possible for supersymmetry to make a 
significant contribution to Rb, if certain conditions on the low-energy supersymmetric 
spectrum are satisfied [B, M. 

In this paper we reexamine this question in the context of a variety of super- 
gravity models, which include the constrained minimal supergravity model considered 
in Ref. ||, as well as non-minimal string-inspired supergravity models. In all cases we 
find R s b nsy ^ 0.0004, which brings the Standard Model prediction at most one-fifth of 
a standard deviation closer to the experimental result. Moreover, R s b usy could well be 
negative, worsening the Standard Model fit, although only slightly. We then turn to 
the minimal supersymmetric standard model (MSSM), where the many parameters 
are a priori independent, and seek to delineate the region of parameter space which 
yields large values of Rb- This exercise confirms the results of Ref. ||, that light 
charginos and top-squarks of definite composition are required. However, we take 
the further step of trying to quantify how large a region of parameter space this is. 
In other words, what kind of fine-tuning is involved in attaining such large values of 
Rb- To this end we apply all known experimental constraints which may affect the 
preferred region of parameter space: (i) that B(b — > sj) be within the allowed CLEO 
range, (ii) that the lightest supersymmetric particle (LSP) be neutral and colorless, 
(iii) that the invisible top-quark branching ratio (i.e., B(t — > iiX? 2 )) be within ex- 
perimental limits, (iv) that the invisible Z width be within LEP limits, and (v) that 
B(Z — > X1X2) be within LEP limits. By sampling a large number of random values 
for the six-plet of parameters that determine R s b usy and the other five experimental ob- 
servables, we conclude that, of all points in parameter space that yield R s b usy > 0.0020 
(la), only a few percent satisfy the combined additional experimental constraints. 
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We should remark that the region of MSSM parameter space where R s b usy is 
enhanced, has recently become the focus of attention for another reason. It has been 
argued that the apparent disagreement between the LEP determination of a s and 
the corresponding value obtained from low-energy measurements, may hint to the 
possibility of new physics M. Moreover, light supersymmetric particles seem to fit 
the bill by shifting downwards the LEP-extracted value of a s [[|, if observables like Rb 
can be brought to better agreement with experiment. Our restrictions on parameter 
space apply also to this case, to the extent that sufficiently large values of Rb are 
called for. 



2 Supersymmetric contributions to 

Besides the one-loop Standard Model contributions to Rb, in supersymmetric models 
there are four new diagrams, as follows: 

• The charged-Higgs-top-quark loop, depends on the charged Higgs mass and 
the t — b — H ± coupling. For a left-handed b quark the coupling is oc m t / tan/? 
whereas for a right-handed b quark it is oc mj, tan (3. Therefore, for smallQ (large) 
tan/3 left- (right)-handed 6-quark production is dominant. (For tan/5 3> 1, the 
value of rrib impacts the contribution significantly.) It has been shown that the 
H ± — t contribution is always negative [M, a fact which makes the prediction 
for Rb in two Higgs-doublet models always in worse agreement with experiment. 

• The chargino-top-squark loop, is the supersymmetric counterpart of the H-t 
loop discussed above. The chargino mass eigenstate is a mixture of (charged) 
Higgsino and wino, and the coupling strength is a complicated matter now 
because it involves the top-squark mixing matrix and the chargino mixing ma- 
trix. However, because only the Higgsino admixture in the chargino eigenstate 
has a Yukawa coupling to the t—b doublet, generally speaking a light chargino 
with a significant Higgsino component, and a light top-squark with a significant 
right-handed component are required for this diagram to make a non-negligible 
contribution to Rb @. 

• The neutralino-bottom-squark loop, is the supersymmetric counterpart of the 
neutral-Higgs-bottom-quark loop. The coupling strength of Xi ~ b ~ b is also 
rather complicated, since it involves the bottom-squark mixing, the neutralino 
mixing, and their masses. However, it can be non-negligible since it is pro- 
portional to mfttan/3 for the left-handed 6-quark. Therefore in the high-tan /3 
region we have to include this contribution. R 



^dn supergravity models, the radiative electroweak breaking mechanism requires tan/? > 1. 
2 This term was not included in our previous study in terms of the parameter p7[ , although it 
was pointed out that its effects were non- negligible for tan j3 3> 1 ■ 
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• The neutral-Higgs-bottom-quark loop, involves the three neutral scalars, two 
CP-even (h,H) and one CP-odd (A). For the h (H ) neutral Higgs boson the 
coupling to the bottom quark is oc m& sin a (cos a) which in the absence of a 
tan/3 enhancement makes its contribution negligible. For the A Higgs boson, 
the coupling to bb is oc tan (3 and the A-dependent contribution can be large 
and positive if tra ^ 90 GeV and tan/3 ^ 30 [|10j. Since we restrict ourselves 
to tan (3 & 20, and uia ^ 100 GeV in supergravity models, this contribution is 
neglected in what follows. 

Our computations of Rl usy have been performed using the expressions given in Ref. |J. 
Even though these formulas are given explicitly, the details are quite complicated 
by the presence of various Passarino-Veltman loop functions. As a check of our 
computations, we have verified numerically that the results are independent of the 
unphysical renormalization scale that appears in the formulas. 



3 Rh in supergravity models 



We consider unified supergravity models with universal soft supersymmetry breaking 
at the unification scale, and radiative electroweak symmetry breaking (enforced using 
the one- loop effective potential) at the weak scale [11]. These constraints reduce the 



number of parameters needed to describe the models to four, which can be taken to 
be m 1( / 2 ,(o = m o/ m i/2, £,a = A/mi/ 2 , tan/3, with a specified value for the top-quark 
mass, which we take to be m t = 175 GeV. Among these four-parameter supersym- 
metric models, we consider generic models with continuous values of m+ and discrete 
choices for the other three parameters: 



tan j3 



2,10,20 ; 



Co 



0,1,2,5; 
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The choices of tan/3 are representative; higher values of tan/3 are likely to yield 
values of B(b — > 57) in conflict with present experimental limits [ I2"| . The choices 
of £0 correspond to Trig « (0.8,0.9, 1.1, 1.9)mg. The choice of A has little impact on 
the results. We also consider the case of no-scale SU(5) x £7(1) supergravity [O. In 



this class of models the supersymmetry breaking parameters are related in a string- 
inspired way. In the two-parameter moduli scenario £0 — £,a = fL3| , whereas in the 

73' ^ A 



dilaton scenario £0 
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A series of experimental constraints and 
predictions for these models have been given in Ref. ||15|| . In particular, the issue of 
precision electroweak tests in this class of models has been addressed in Refs. |16 , [17]. 



The predictions for i?™ sy in the generic supergravity models are shown in 
Fig. [I]. Only curves for tan /3 = 2, 10 are shown, since the corresponding sets of 
curves for other values of tan /3 fall between these two sets of curves. The results in 
SU(5) x £7(1) supergravity, for the same values of the parameters, differ little from 
those in the generic models. To add generality to our result and to compare with 
the study made below in the case of the MSSM, we have also considered the case 
where the four supergravity parameters are chosen at random. We have sampled 
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10,000 random four-plets of parameter values, in the ranges: 50 < m\/2 < 500 GeV, 
< £ < 5, -5 < U < 5, and 1 < tan/3 < 40. As expected, we find R s b nSY < 0.0004. 
A histogram of the relative distribution of Ft b usy values is shown in Fig. |2], which shows 
that R s b nsj is equally likely to be positive or negative, and that the "preferred value" 
is very small. 

In almost all cases the largest positive contribution to i?™ sy comes from the 
chargino-top-squark loop. As expected, the largest contribution from this diagram 
happens for points with the lightest chargino masses (which correspond to the lightest 
t\ masses) since supersymmetry is a decoupling theory. However, even the largest 
value (~ 0.0004) is still very small compared with the largest possible result in a 
generic low-energy supersymmetric model. The reason for this is that while the 
smallest possible chargino and a top-squark masses are required for an enhanced 
contribution, it is also necessary that the top-squark be mostly right-handed and that 
the chargino has a significant Higgsino component |J. A scatter plot of the values of 
these couplings is shown in Fig. |3|, where the right-handed component of the lightest 
top-squark (T 12 ) is plotted against the higgsino component of the lightest chargino 
(V12). The first requirement is in fact attainable in these models {i.e., \Tx%\ ~ 1), but 
the former is not {i.e., | V12I ~ 0.4). This behavior is expected in supergravity models 
with radiative electroweak symmetry breaking since for light charginos \fi\ 3> M2, 
which makes the lightest chargino mainly a wino instead of a Higgsino. (This is also 
the reason why the results in SU{5) x U{1) supergravity differ little from those in 
the generic models.) Concerning the other contributions to R s b usy , the charged-Higgs- 
top-quark loop is always negative, and is enhanced for either small or large values of 
tan /3. The neutralino-bottom-squark contribution is almost always smaller than the 
chargino-top-squark contribution and not of definite sign. 

Thus, Rb in supergravity models could improve the Standard Model fit to the 
LEP data by at most one-fifth of a standard deviation, and it could well worsen the 
fit by the same small amount. 

4 R b in the MSSM 

We now relax the supergravity assumptions that correlate the various supersymmet- 
ric parameters, in order to explore the region of the MSSM parameter space that 
yields large values of Rb. This region is then subjected to all available experimental 
constraints, which have the effect of restricting it significantly. The R s b usy observ- 
able depends on six basic parameters: the elements of the chargino mass matrix 
(M 2 , fi, tan/3) which determine the chargino masses and their composition, and the 
two top-squark {ti,t 2 ) masses and their mixing angle (# t ~). In addition, the Higgs 
boson masses enter, although for tan j3 < 30 the neutral Higgs boson contribution is 
not relevant ||. The charged Higgs mass is relevant for small values of tan/3, but this 
contribution to Rb is always negative and will be neglected in what follows {i.e., we 
take a large charged Higgs mass). Neglecting altogether the Higgs-boson contribution 
to R b wy is generally justified when looking for the largest values of Rb, although some 
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exceptions exist for rather low values of the pseudoscalar Higgs boson mass (m^) and 
large values of tan j3 J7| . 

We have sampled a large number of random choices for the six-plet of param- 
eters, in the ranges < M 2 , m^m^ < 250 GeV, 1 < tan/? < 5, and < 9i < ir. 
We have concentrated on the small-tan f3 region since for large values of tan/5, where 
R s b nsy may also be enhanced, B(b — > 57) is likely to exceed significantly the allowed 



CLEO range ||12|| . A histogram depicting the relative distribution of calculated values 
of R s b usy is shown in Fig. f|. [Note that small negative values of R s b usy are possible 
(although not very likely).] As expected, there is a steady decline in the likelihood of 
the larger values of R s b nsy , with the following relative distributions 

Condition Comment /i > /i < 



R s b usy > 0.0008 Ok at 95%CL 12% 



O.O70 



R s b usy > 0.0014 Ok at 90%CL 3.0% 1.9% (2) 
R s b usy > 0.0020 Rf M up by la 0.87% 0.52% 
R s b nsy > 0.0027 Ok at la 0.20% 0.12% 

where for example, only ~ (2 — 3)% of the sampled six-plets yield Rl nsy > 0.0014, 
which brings the total Rb prediction inside the 90%CL allowed range for m t = 
175 GeV. Also, to be within the la range, and thus to significantly improve the 
fit to the LEP data, requires considerable fine tuning. If the sampling region were to 
be extended (i.e., a larger mass interval or larger values of tan/3) we would find that 
the new points fall into the lowest bins in Fig. [|, while the population of the bins with 
large Rb values would correspondingly decrease. The distributions in Eq. (g) would 
change accordingly. 

In Fig. [5] we show the maximum attainable value of R s b usy as a function of 
the chargino mass, which makes apparent the need for a light chargino if R s b usy is to 
be enhanced. In fact, the maximum value is obtained for m ± — > \Mz (and also 

— > ^Mz). This phenomenon has been observed before in the context of the e& 
approach to the problem [H |i~6| , and corresponds to a wavefunction renormalization 
effect when the particles in the loop go on shell. 

For concreteness, in what follows we concentrate on the ^ 1% of the points 
that increase Rb by at least la; the other cases in Eq. (|) depend on the choice of 
m t . By design, this sample contains 1000 points for each sign of fi. First let us show 
that the region of parameter space which leads to enhanced values of R s b usy is indeed 
characterized by light higgsino-like charginos and light right-handed-like top-squarks. 
In Fig. H we show the distribution of points in the (m x ± , myj plane that lead to i?™ hy > 
0.0020, and in Fig. [7] we show the corresponding distribution in (|Vi 2 |, |T\ 2 1 ) space. If 
R s b usy > 0.0020 is indeed required, then LEPII should see the lightest chargino and 
possibly also the ligthest top-squark. The observed marked preference for large values 
of the V12, T12 admixtures becomes evident when one considers the top-quark Yukawa 
coupling 

A* QzHt R — > A 4 bintR , Ai tLntR , (3) 



5 



Table 1: The 95%CL upper limit on the invisible top-quark branching ratio, obtained 
from a comparison of DO data with theoretical expectations. All masses in GeV, all 
cross sections in pb. 



m t 


<X D0 
mm 


th 

max 


min 


Dt— >inv 
max 


150 


3.6 


13.8 


0.51 


0.49 


160 


3.4 


9.53 


0.60 


0.40 


170 


3.1 


6.68 


0.68 


0.32 


180 


2.5 


4.78 


0.72 


0.28 


190 


2.2 


3.44 


0.80 


0.20 


200 


2.1 


2.52 


0.91 


0.09 



which picks out ti — > Ir, xt ~^ H ± , \\ ~^ H°. This interaction also leads to an 
enhanced coupling between the top-quark, a higgsino-like neutralino, and a right- 
handed top-squark, and may lead to enhanced exotic decays of the top-quark, as we 
discuss below. Comparing Fig. |3| with Fig. |7| we see why R s ^ sy is always suppressed in 
supergravity models: the regions in (IV12I, |2~i2|) space are completely non-overlapping. 
We should note that for // > it is easier to obtain large-i?™ sy solutions since the 
chargino can have a larger higgsino component (as Fig. [7] shows), i.e., [ V12I — I cos^> + | 
with tan(20 + ) oc — /x cos (3 + M 2 sin j3 [Tj| • 

The parameters that determine R s ^ sy also determine the following observables: 

1. B(b — > s-y), as measured by CLEO to be (1 — 4) x 10~ 4 fll9f , has been computed 



as in Ref. JT^|. The charged-Higgs loop is negligible in the limit of a large 
charged Higgs mass (as assumed above) and in general its contribution is not 
large and does not significantly affect whether B(b — > &y) is in agreement with 
experiment or not. 

2. m x o < {m x ±,mi^\, is the cosmological requirement of a neutral and colorless 
lightest supersymmetric particle (LSP). The lightest neutralino mass follows 
from the inputs to the chargino mass matrix and the usual assumption relat- 
ing Mi to M.2- This constraint has already been included in Eq. (Q) and in 

Figs, am 

3. Bit — > £1X12) leads to an "invisible" decay of the top quark if the t\ bxf — * 
b(qq', £vi)xi decay products do not pass the standard top-quark cuts, as may 
likely be the case given the larger amount of missing energy and the softer 
leptons that accompany this decay. Taking the 95% CL lower bound on the 
top-quark cross section from DO (see Fig. 3 in Ref. [p0|1), and dividing it by 
the "upper" estimate of the theoretical cross section (see Table 1 in Ref. [pT|), 
one can obtain a 95%CL lower bound on [B(t — > VT6)] 2 . As a function of 
mt, this exercise is carried out in Table [I]. For m t = 175 ± 15 GeV we obtain 

cr = o.3 T o.i. 
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Table 2: The fraction of MSSM parameter space with R s b usy > 0.0020, where the five 
experimental constraints are satisfied individually, and the combined allowed region, 
for B^iT = 0.3 ± 0.1. 



n > 



n < 



Dt— >inv 
max 




LSP 


^t— >inv 


"pinv 
1 Z 


B Z ~*M 


All 


0.4 


20% 


98% 


41% 


88% 


72% 


5.3% 


0.3 


20% 


98% 


19% 


88% 


72% 


1.7% 


0.2 


20% 


98% 


8.0% 


88% 


72% 


0.9% 


Dt— »inv 
max 


gb^s~f 


LSP 


^t— >inv 


pinv 
1 Z 




All 


0.4 


37% 


97% 


31% 


98% 


95% 


8.9% 


0.3 


37% 


97% 


14% 


98% 


95% 


1.9% 


0.2 


37% 


97% 


8.4% 


98% 


95% 


0.6% 



4. The invisible Z width Y[Z — > XiXi) follows from the neutralino mass and 
composition. It should not exceed 7.6 MeV . 



5. The branching ratio BiZ — > X1X2) a ^ so follows from the neutralino mass matrix, 
and should not exceed 10~ 4 |23l. 



We now impose these additional constraints on the points in MSSM parameter 
space with R s b Uby > 0.0020. First let us give the fraction of parameter space where 
these constraints are satisfied individually, as well as the combined allowed region. 
In Table ^| we carry out this exercise for 5 t-nilv = 0.3 ±0.1, which corresponds 
to m t = 175 =F 15 GeV. One can see that only a few percent of the points with 
interestingly large values of R s b usy are experimentally allowed. Table |2] also shows that 
the b — > 57 and _g t->inv constraints are the most restricting ones, and that they are 
both satisfied only in a small region of parameter space. This trend is shown more 
concretely in Fig. ||, where the distribution in [_B t_+mv , B(b — > 57)] space is shown, 
together with the present experimental limits on these observables. In this figure, for 
ji > (/x < 0): ~ 20 (37)% of the points fall inside the allowed B{b — > 57) region, 
~ 41 (31)% of the points fall inside the allowed B^' mv area, and ~ 5 (9)% fall in 
the combined allowed area. As remarked above, in regions of parameter space with 
enhanced values of , we may also expect enhanced values of B 1 ^™. In Fig. [5] 
we show the distribution of i?^ usy versus B^ irw calculated values (with a cutoff of 
R s b Uby > 0.0020), which indeed shows a distinct correlated enhancement of these two 
observables. 

Finally, we consider the case where all experimental constraints are applied 
simultaneously. The resulting allowed points in the (m x ±, rajj plane, for B^™ = 0.4, 
are shown in Fig. TU. This figure is to be contrasted with Fig. Bl where none of 



the experimental constraints were applied (except for the rather unrestricting LSP 
constraint). We see that the extent of the allowed area is reduced: m™± x : (70 
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60) GeV and m? iax : (110 — > 80) GeV. This restricted parameter space should be 
even more straightforwardly tested at LEPII, with now guaranteed sensitivity to the 
ligthest top-squark. If we strengthen the 5 t_>mv constraint to B^^ v = 0.2, the 
more relevant shrinking of parameter space occurs in the and tan/3 directions: 
m t ~ 2 < 100 GeV and tan/3 ^1.5 are now required. Such strong restrictions would also 
allow detection of the lightest Higgs boson at LEPII, and both top-squarks at the 
Tevatron and LEPII. 



5 Conclusions 

We have studied the contributions to Rb that are expected in low-energy supersym- 
metric models. In the case of supergravity models with radiative electroweak symme- 
try breaking, R s b nsy could improve the Standard Model fit to the LEP data by only a 
small fraction of a standard deviation, and it could well worsen the fit by the same 
small amount. In the MSSM, light higgsino-like charginos and light right-handed-like 
top-squarks are required to obtain a sizeable R s b nsy value. Such values occur in a small 
fraction of the parameter space, which is largely accessible to LEPII searches. This 
region of parameter space is further constrained by five experimental observables, 
most importantly B(b — > 57) and the invisible top-quark branching fraction. Impos- 
ing these additional constraints reduces the allowed parameter space to a few percent 
of its unconstrained size, and makes experimental exploration of this scenario in both 
chargino and top-squarks assured at LEPII. 
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Figure 1: The supersymmetric contribution to Rb as a function of the chargino mass 
in generic supergravity models with tan/3 = 2, 10, £ — — 5, and A = 0. Curves for 
other values of tan (3 fall between the two sets shown. 
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Figure 2: Relative distribution of calculated values of R s b nsy in a sample of 10,000 
random choices for the four-plet of supergravity parameters. 
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Figure 3: The right-handed component of the lightest top-squark (T 12 ) versus the 
higgsino component of the lightest chargino (V12) in a random sample of supergravity 
models. Large values of both T i2 and V\i are required for an enhanced value of it^ usy . 



13 




Figure 4: Relative distribution of calculated values of R s b nsy in a large sample of 
random choices for the six-plet of relevant MSSM parameters. For yU>0 (/x<0), only 
« 0.9% (0.5%) of the points yield R s ™ y > 0.0020 (la). 
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Figure 5: Maximum attainable value of Ft b nsy in the MSSM as a function of the 
lightest chargino mass. All other variables have been integrated out. The broken line 
is a figment of the finite sample size. 
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Figure 6: The distribution of lightest chargino and top-squark masses for points in 
MSSM parameter space with R™ sy > 0.0020 (la). 
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Figure 7: The right-handed component of the lightest top-squark (T 12 ) versus the 
higgsino component of the lightest chargino (V12) for the fraction of MSSM sampled 
points with R™ sy > 0.0020 (la). 
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Figure 8: The calculated value of B(b — > 57) versus ,B t_>mv for points in the MSSM 
parameter space with i?^ usy > 0.0020. The arrows point into the experimentally 
allowed region. For /i>0 (//<0) only 5% (9%) of the points fall inside the combined 
allowed region. 



18 



/z>0 



Rl US7 > 0.0020 



/x<0 



0.005 i i i i i i i i i i i i i i i i i i i i i i i i i 



0.004 — 



j^susy 



0.003 — 



0.002 



0.001 















0.2 



0.4 0.6 
gt-»inv 



0.8 



i i i i l i i i i i i i i i i i i i i i 



••■■'•-^'. ; "; ; -vV"- 











0.2 0.4 0.6 



0.8 



B 



t-»inv 



Figure 9: The calculated value of R s b usy versus £ Winv for points in the MSSM pa- 
rameter space with R s ^ sy > 0.0020. Note the correlated enhancement of these two 
observables. 
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Figure 10: The distribution of lightest chargino and top-squark masses for points in 
MSSM parameter space with i?™ sy > 0.0020 (lcr), when all experimental constraints 
are imposed (with B^£™ = 0.4). Because of the finite sample size, the whole region 
spanned by the shown points should be considered viable. 
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